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ABSTRACT 

We carry out an extended analytic study of how the tilt and faster-than-radial 
expansion from a magnetic field affect the mass flux and flow speed of a line- 
driven stellar wind. A key motivation is to reconcile results of numerical MHD 
simulations with previous analyses that had predicted non-spherical expansion 
would lead to a strong speed enhancement. By including finite-disk correction 
effects, a dynamically more consistent form for the non-spherical expansion, and 
a moderate value of the line-driving power index a, we infer more modest speed 
enhancements that are in good quantitative agreement with MHD simulations, 
and also are more consistent with observational results. Our analysis also explains 
simulation results that show the latitudinal variation of the surface mass flux 
scales with the square of the cosine of the local tilt angle between the magnetic 
field and the radial direction. Finally, we present a perturbation analysis of the 
effects of a finite gas pressure on the wind mass loss rate and flow speed in both 
spherical and magnetic wind models, showing that these scale with the ratio 
of the sound speed to surface escape speed, a/v esc , and are typically 10-20% 
compared to an idealized, zero-gas-pressure model. 

Subject headings: MHD — Stars:winds — Stars: magnetic fields — Stars: early- 
type — Stars: mass loss 

1. Introduction 

In a recent paper (ud-Doula & Owocki 2002; hereafter UO-02), we presented numerical 
magnetohydrodynamic (MHD) simulations of the effect of a stellar dipole magnetic field 
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on the line-driven stellar wind from a non-rotating, hot, luminous star. We showed that 
the overall influence of the field on the wind depends largely on a single, dimensionless, 
'wind magnetic confinement parameter', r]* (= B^ q Rl/ Mv^), which characterizes the ratio 
between magnetic field energy density and wind kinetic energy density. Because the field 
energy declines faster than the wind energy, in the outer regions the magnetic field is always 
dominated by the radial wind outflow, which thus asymptotically stretches the field into a 
radial configuration, regardless of the strength of 77*. For weak confinement 77* < 1, this 
radial opening of the field applies throughout the whole computational domain. But for 
stronger confinement 77* > 1, the magnetic field remains closed near the surface and over a 
limited range of latitude around the magnetic equator. 

In this paper we provide further analysis and interpretation of how key wind properties, 
namely the surface mass flux and asymptotic flow speed (see figs. la,b in §2 below), are 
modified by the presence of a magnetic field. One central motivation is to reconcile results of 
our numerical MHD simulations with the previous scaling analysis done by MacGregor (1988; 
hereafter M-88), which had predicted that the faster-than-radial divergence of magnetic flux 
tubes would lead to substantial (factor 3-4) increase in the terminal velocity of a line-driven 
wind, compared to that of the non-magnetic, spherical-expansion case. Since the implied flow 
speeds ranging up to 5000-6000 km/s are never observed through blue edges of UV P Cygni 
lines from hot stars (Prinja & Howarth 1986; Howarth and Prinja 1987), that prediction 
could be a basis for inferring that hot stars lack magnetic fields with sufficient strength to 
substantially influence their stellar winds. 

Here we develop (§3) a simplified formulation of the one-dimensional (ID) equations for 
a steady line-driven wind. Applying this to cases with faster-than-radial area expansion (§4), 
we show that the prediction of a 3-fold increase in flow speed is the consequence of certain 
assumptions and idealizations (point-star radiation, phenomenological flow expansion, and 
large line power-index a) in the previous M-88 analysis; when these assumptions are relaxed 
(to include finite-disk correction, a dynamical flow divergence, and more realistic a) the 
expected effect on flow speeds is consistent with the more modest ^ 50% increase typically 
found in our MHD simulations. We further show that the inferred simulation scaling of the 
base mass flux with the square of the cosine of the surface field tilt angle can be understood 
from a simple modified ID tilted flow analysis (§5). After a brief discussion of the implications 
for stellar wind structure (§6), we conclude with a summary (§7) of our main results. 
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2. Summary of Key Results of 2D MHD Simulations 

Figures la,b summarize two key results of our 2D MHD simulations, wherein we study 
the dynamical competition between magnetic field and wind by inserting a dipole magnetic 
field of varied strength (parameterized by 77*) into a previously relaxed spherically symmetric 
CAK wind. 1 For the asymptotic states when the regions of magnetic field opened by the 
wind have a nearly steady-state flow, both figures show results for three cases, with the 
degree of magnetic confinement ranging from weak (77* = 0.1), through moderate (77* = 1), 
to strong (77* = 10). 

In fig. la, the dotted curves show the latitudinal variation of the base mass flux, while 
the solid curves show this same flux scaled by where \ib is the radial projection of a unit 
vector along the base magnetic field. The near constancy of this scaled flux is explained in 
§5, which uses a simple tilted-flow analysis of the base line-driving to show that the radial 
mass flux does indeed vary with the square of the field's radial projection cosine, M ~ \j? B . 

At the magnetic poles, the radial orientation of the field means the flow is also radial, 
with essentially the same base mass flux as in the spherical model. 2 

Fig. lb plots the latitudinal variation of the flow speed at the maximum model radius 
Rmax = 6R*, scaled by the corresponding speed in the spherical, non-magnetic (B = 0) 
model. The lower speeds toward the magnetic equator (fx « 0) reflect the various rela- 
tive degrees of equatorial flow compression, which increases the density, and so reduces the 
radiative acceleration. In the highest magnetic confinement case 77* = 10, the increased 
speeds at mid-latitudes (0.1 < < 0.8) result from a combination of the reduced mass flux 
and greater flow expansion, both of which reduce the density, and so increase the radiative 
acceleration. 

While such latitudinal trends thus make sense qualitatively, our focus below will be on 
developing a semi- quantitative understanding of the wind properties at the magnetic poles, 
where the radial orientation of the flow makes possible a ID analysis, if one takes proper 
account of the non-spherical expansion of the flow area. This was essentially the approach 



1 The simulation results given here are analogous to those presented by UO-02 and ud-Doula (2002), 
but correct a minor coding error that lead the polar flow to have an artifically enhanced mass flux and an 
artificially reduced flow speed. As in UO-02, the model here is assumed isothermal, with a constant sound 
speed a = 26.3 km/s, implying a ratio s = a 2 jv\ sc = 0.0014 between the gas internal energy and surface 
escape energy. (See §3 and the Appendices.) 

2 The few percent increase is a second-order effect, associated with the finite sound speed and the resulting 
detachment of the flow critical point from the surface, where faster areal expansion reduces the density, and 
so allows driving of a somewhat greater mass flux. See Appendix B. 



-4- 




Fig. 1. — a. The latitudinal variation of the base mass flux in open field regions of 2D MHD 
simulations. The dotted curves represent the mass flux scaled by the corresponding spherical 
model. The solid curves add an additional scaling by ji 2 B , where hb = -£>* • f is the radial 
projection of a unit vector along the base magnetic field, b. Latitudinal variation of the flow 
speed at the maximum radius R m ax = QR* of these same 2D MHD simulations, scaled by 
the corresponding speed in the spherical, non-magnetic (B = 0) model. 
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taken by M-88, but the factor 3 or more speed enhancements found there are substantially 
greater than those in fig. lb, which even for the strongest magnetic confinement r]* = 10 
shows a polar speed enhancement of only 1.54. The next two sections develop and apply a 
formalism for understanding this difference. 



3. Flow Analyses for Spherical Expansion 

3.1. CAK Model of a Wind Driven by a Point-Source of Radiation 

To provide a basis for our analysis of how the tilt and areal divergence of a magnetic field 
modify a wind outflow, let us first recast the basic Castor, Abbott and Klein (1975; hereafter 
CAK) formalism in a somewhat more tractable, physically transparent form. First, for a 
steady, !D stellar wind driven radially by line-scattering of radiation from a central point- 
source, the associated equation of motion for the radial velocity v as a function of radius r 
takes the form 



2 

f / 

1 - 



dv GM{\ - T) 2a 



v' 2 



2 



V-r- = j + 9cAK + ' W 

dr r A r 



where GM(1 — T)/r 2 is the effective gravitational acceleration, and the Eddington parameter 
T = K, e L/47rGMc accounts for the effect of electron scattering opacity K e interacting with 
the stellar luminosity L, with c the speed oflight. The square-bracket factor with the sound 
speed a on the left-hand-side accounts for gas pressure effects that allow for smooth mapping 
of the wind outflow onto a hydrostatic atmosphere below the sonic point, where v — a. On 
the right-hand-side the flow's areal expansion also gives rise to a pressure-form term, 2a 2 /r, 
which compared to the competing gravitational term is of order s = (a/v esc ) 2 ~ 0.001, 
where v esc = a/2GM(1 — T)/R* is the effective escape speed from the stellar surface radius 
R*. The outward driving against gravity through the sonic point and beyond must thus 
be principally via the line acceleration, which under CAK's original assumption of radially 
streaming radiation from a point source can be written as 

1 n e LQ ( dv/dr\ a 
9 cak = jz -r -r—T- —7s— (2) 



a) 4Trr 2 c \pcQn, 

where p is the density, and a is the CAK exponent. Here the line normalization Q is related to 
the usual CAK "k" parameter through k = Q x ~~ a (v th /c) a /(l — ck), but offers the advantages 
of being a dimensionless measure of line opacity that is independent of the assumed ion 
thermal speed v t h, with a nearly constant characteristic value of order Q ~ 10 3 for a wide 
range of ionization conditions (Gayley 1995). 
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To proceed let us define the gravitationally scaled inertial acceleration 

r 2 v dv/dr 



w 



GM(1 - T) ' 



(3) 



In terms of an inverse radius coordinate x = 1 — R*/r, note that w' = dw/dx, where 
w = v 2 /v 2 sc represents the ratio of wind kinetic energy to the effective gravitational binding 
v 2 sc /2 = GM(1 — Y)/R from the stellar surface radius R*. We can then rewrite the equation 
of motion (1) in the dimensionless form 



(1 - s/w) w' = -l + C(w') a + 2s/(l - x) 



(4) 



where s = a 2 /v 2 sc , and we have eliminated the density p in favor of the mass loss rate 
M = Anr 2 pv for the assumed steady, spherical expansion, with the line-force constant defined 
by 



C = 



a 



Mc 2 



QT 



(5) 



Note that, for fixed sets of parameters for the star (L, M, T) and line-opacity (a, Q), this 
constant scales with the mass loss rate as C ~ 1/M a . 

The essential elements of interest in this paper can be analyzed without including the 
small gas pressure terms, and so we shall henceforth consider only the limit of vanishing 
sound speed a ~ y/s — > 0. (The corrections from a small but finite s are examined in the 
Appendices.) For this limit, fig. 2 illustrates a graphical solution to eqn. (4) for various 
values of C. For high M or small C there are no solutions, while for small M or high C 
there are two solutions. The CAK critical solution corresponds to a maximal mass loss rate, 
and requires a tangential intersection between line-force and combined inertia plus gravity, 
for which 

1 , (6) 



a C c w' c 



a-l 



and thus 



with 



a 



1 — a ' 
1 



Q — 

c a a {l-a) 1 - a 



(7) 



(8) 



Applying this in eqn. (5), we can then solve for the mass loss rate, yielding the standard 
CAK scaling 



M, 



a 



CAK 



a 



' QT ' 

i-r 



(l-a)/a 



(9) 
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Fig. 2. — Graphical solution of the dimensionless equation of motion (4) representing a ID, 
point-star, zero-sound-speed CAK wind, as controlled by the constant C ~ 1/M a , with 
determine the intersections of the line-force Cw' a (solid lines) with the combined gravity 
plus inertia 1 + w' (dot-dashed lines). If M is too big, there are no solutions; if M is small 
there are two solutions. A maximal value M = Mcak defines a single, "critical" solution. 
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Moreover, since the scaled equation of motion (4) has no explicit spatial dependence, the 
scaled critical acceleration w' c applies throughout the wind. This can thus be trivially inte- 
grated to yield w = ax /{I — a), which is equivalent to the CAK velocity law 

v(r) = Voo \ K l-^j , (10) 

where the wind terminal speed = v esc ^J a / (I — a) scales with v esc = \J2GM{1 — T)/R*, 
the effective escape speed from the stellar surface radius R*. 



3.2. The Finite-Disk Correction Factor 

The above analysis has so far been based on the idealization of radially streaming radi- 
ation, as if the star were a point-source at the origin. This was the basis of the original CAK 
wind solutions, although they did already identify (but did not implement) the appropriate 
"finite-disk correction factor" to account for the full angular extent of the star [See CAK 
eqn. (50).], 

f _ (1 + *) 1+a ~ (1 + ^l) 1+a (u) 
Jd (l + a)a(l + (7)°(l-^) ' 1 } 

where a = (r/v)dv/dr — 1 and = 1 — Rl/r 2 . When this factor is included to modify the 
point-star CAK line-acceleration (2), its complex dependence on radius, velocity, and velocity 
gradient complicates somewhat the solution of the full equation of motion (1). But full 
solutions have been derived independently by Friend and Abbott (1986) and Pauldrach, Puis, 
and Kudritzki (1986). They yield a somewhat reduced mass loss rate Mfd ~ Mcak/{^ + 
a) 1 /", higher terminal speed ~ 3t> esc , and flatter velocity law, roughly characterized by 
replacing the exponent 1/2 in eqn. (10) with a somewhat higher value (5 ~ 0.8. 

If one simply assumes such a canonical "beta" velocity law, fig. (3) illustrates the specific 
spatial variations of fd{x) for various CAK exponents a and velocity exponents f3. All these 
varied cases show a quite common general behavior. For the initial acceleration from a static 
stellar surface (r « or x — > 0; where v — > with finite dv/dr, implying a — > oo), fa has a 
starting value fa* ~ 1/(1 + a). It then increases outward, past unity at the point where the 
wind has a locally isotropic expansion (with dv/dr = v/r, so that a — 0), and then finally 
falls back asymptotically toward unity at large radii (x — > 1), where the star is indeed well 
approximated by a point source. 

To understand the dynamical effect of this finite-disk factor, and indeed to see how it 
interacts with effects of a magnetic field, it is helpful to examine how such a typical radial 
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Fig. 3. — Spatial variation of the finite-disk-correction factor fd, plotted vs. scaled inverse 
radius x = 1 — R*/r, for various CAK exponents a and velocity-law exponents j3. The heavy 
solid curve shows our standard case with a = 0.6 and (3 = 3/4. The dashed curves use 
(5 = 1/2 and the dot-dashed curves (5 — 1, and are each plotted for CAK exponents a = 0.5 
and a = 0.7 that span the nominal range. A key point is that all these varied curves have a 
similar form, rising sharply from a base value [fa* = 1/(1 + a)], crossing above unity (dotted 
line), and then relaxing back to it at large radial distances, x — > 1. 
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variation influences the scaled acceleration w' through solutions to the appropriately modified 
equation of motion [cf. eqn (4)], 

w' = -l+fC*(w') a , (12) 

where for convenience we have defined scaled quantities / = fd/fd* and C* = fd*C, and, 
for simplicity, again taken the zero-sound-speed limit s — 0. Since / increases outward 
from unity at the surface, the surface radius r = R* (x — 0) now represents a critical point 
(also known as a choke or throat) that fixes the required minimal value of = C c = 
l/a a (l — a) 1-0 to allow an accelerating solution from this wind base. Accounting then for 
the mass loss rate scaling of ~ l/M", we readily see that the net result is to reduce the 
maximal allowed mass loss to 

M fd = fdl a M C AK = M CAK /(1 + af a . (13) 

This thus provides a simple rationale for this aspect of the detailed numerical solutions of 
Friend and Abbott (1986) and Pauldrach et al. (1986). 

As the finite-disk factor increases from its reduced base value, the lower mass load- 
ing allows for stronger CAK line-force, with now the resulting greater acceleration being 
"leveraged" by its effect in further increasing the line-force. From the graphical solution 
in fig. 2, one sees that a modest increase in the line-acceleration term [C(it/) a ] above the 
critical solution (with C = C c ) leads to a solution with a much larger scaled acceleration 
w'. [There also arises a much shallower acceleration, but as discussed by Feldmeier et al. 
(2002), maintaining this beyond the critical point requires a special kind of outer boundary 
condition, whereas the more rapid acceleration is compatible with the standard assumption 
of simple expansion into a vacuum.] 



Fig. 4 illustrates this strong sensitivity to the finite-disk correction in terms of vto', 
which sets the scaling of terminal speed Voo ~ v esc \fw' . Note how a modest, less-than-factor- 
two change in / leads to large increases in this speed scale \fw', with the trends being most 
dramatic for the largest a. From eqn. (12), we see that in the limit of large acceleration 
w' >> 1, the speed scaling approaches the form 

v^7~ (C c f) 1/2(1 - a) ; (14) 

which is illustrated through the dashed curves in fig. 4. For the specific case a — 1/2, we 
can use the quadratic formula to obtain a simple expression that applies to the full range, 

vV=/[l + v / l-l// 2 ] , (15) 
where the choice of the "plus" root gives the steeper of the two acceleration solutions. 
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Fig. 4. — Root of the scaled acceleration y/v? that solves the equation of motion (12), plotted 
as a function of the surface- normalized finite-disk factor /. Since the wind terminal speed 
scales as Voo ~ \fw'v esc , note that even modest increases in the finite-disk factor can lead to 
a very high speed wind (especially for larger a). The dashed curves denote the asymptotic 
solutions ~ (C c f) l ^ 2t < l ~ a ^ that apply in the limit of large acceleration w' ^> 1. 
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Fig. 5. — Comparison of the full finite-disk solutions (solid) with integrations that use a beta- 
velocity law to obtain spatially explicit approximations for the finite-disk factor (dashed). 
Left: The ratio of terminal to escape speed vs. CAK exponent a. Right: Radial variation 
of the velocity for the a = 0.6 case. In both panels, results for a model using (3 = 3/4 
to evaluate the finite-disk factor are also included, but are indistinguishable from the full 
solution plotted by the solid line. The dotted curves show corresponding results for the 
point-star CAK model. 
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To obtain specific values for the flow speed, we can numerically integrate the scaled 
equation of motion (12) accounting for the complex functional dependence of the finite disk 
factor (11) on radius, velocity, and velocity gradient. The left panel in fig. 5 shows how the 
ratio of terminal to escape speed depends on the CAK power index a; the right panel plots 
the full radial variation of velocity for the a = 0.6 case. The dashed curves compare results 
for a similar integration when a canonical beta velocity law is used to evaluate the finite-disk 
factor. For (3 = 0.75 the curves are indistinguishable from the full calculation, and other 
values of (5 also give very good approximations. 

A basic conclusion is thus that one may derive accurate finite-disk-corrected wind solu- 
tions using a beta velocity law to obtain spatially explicit approximations of the finite-disk 
factor, with especially good agreement obtained for (5 ~ 0.75. 



4. The Effect of Non-Spherical Area Expansion 

4.1. Kopp-Holzer Expansion in a Point-Star CAK Wind 

The above formulation of standard line-driven wind theory provides a convenient basis 
for understanding the results of models extended to include magnetic fields. As first analyzed 
by MacGregor (M-88), one potentially key effect regards the faster-than-radial expansion of 
the area in a flow tied to a rapidly diverging magnetic field. Following Kopp & Holzer (1976), 
it is convenient to consider an infinitesimal flow tube with a non-radial expansion, but still 
with a radial orientation, such as would apply, for example, to flow emanating from the poles 
of a surface magnetic dipole. The non-radial expansion factor can then be written as 

, , . Air) Rl B(R*) Rl , x 

hir) = -i- 1 -^ = -^r^f, 16 
v ; A, r 2 B(r) r 2 ' v ; 

where the latter equality follows from the V • B = property of any magnetic field, since 
that then implies d(BA)/dr = and so A ~ l/B along such a radial flow tube. 

In the absence of a detailed dynamical solution for the magnetic field divergence, M-88 
adopted the Kopp & Holzer (1976) phenomenological non-spherical expansion factor 3 , 

, _ h max exp[(r - Ri)/a] + 1 - (h max - 1) exp[(iZ* - Ri)/a] 

exp[(r - Ri)/o-\ + 1 



3 We use the notation h instead of Kopp and Holzer's / to distinguish from the finite-disk correction 
defined above. 
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where Ri, a, and h max are free parameters, the first two specifying the location and extent of 
the non-spherical expansion region, and the last fixing the total overall level of non-spherical 
expansion. Following M-88, we fix sample values of the two Kopp-Holzer parameters, R± = 
1.25-R* and a — 0.1, and compare models with various h max = 1.5 and 2. 

Applying this Kopp & Holzer expansion factor to CAK line-driven wind models com- 
puted in the idealized point-star approximation, M-88 found substantial (more than factor 
three) increases in wind terminal speed for even modest values of the total non-spherical 
expansion, h max ~ 2. The above formalism provides a convenient way to understand those 
results, and indeed to extend them to account for the combined effects of a non-spherical 
expansion and a finite-disk correction. (See §4.2.) Noting from eqn. (2) that the line ac- 
celeration scales with the density p as Qcak ~ 1/P a > we see from the steady-state mass 
conservation relation p ~ 1/vA ~ h(r)/vr 2 that along a radial flow with non-spherical 
expansion ( "nse" ) the line force has the form 

9nse = h(r) a g CA K • (18) 

The relevant equation of motion thus takes the same general form as given in eqn. (12) if 
we simply replace the finite-disk correction with the non-spherical expansion, / — > h a . With 
this redefinition, the dependence of the scaled acceleration w' on h a is thus again given by 
fig. 4. 

Fig. 6a illustrates results for this case, showing the spatial variation (in terms of the 
inverse radius coordinate x) of logw' (left) and wind speed ratio y/w = v/v esc (right). (The 
latter is obtained by straightforward numerical integration of the acceleration w' over x.) In 
both panels, the dotted and dashed curves represent respectively the point-star and finite- 
disk solutions for spherical expansion, while the dot-dash and solid curves respresent the 
corresponding non-spherical solutions for this standard M-88 case with a = 0.7 and h max = 2. 
Note that both the acceleration and velocity laws in this point-star, non-spherical expansion 
case (dot-dash curve) agree roughly to what is found in the corresponding spherical, finite- 
disk model (dashed curve). In this sense, the M-88 results for this h max can be roughly 
thought of as simply substituting the non-spherical expansion for the finite-disk correction. 

For this case with h max = 2 we find in particular a terminal speed ratio Voo/v esc = 
5.2, which, for the M-88 assumed stellar parameters that set the effective escape speed 
v esc = 974 km/s, yields a terminal speed « 5145 km/s, or slightly above the value 
= 4799 km/s quoted in table 1 of M-88. (The modest difference is due to the neglect 
here gas pressure terms that were included in M-88's more elaborate analysis; see Appendix 
A and B.) When compared to the terminal speed expected for point-star spherical expansion 
foo = a/«/(1 — ct)v esc = 1.53 x 974 = 1488 (which is slightly below the value 1512 km/s 
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Fig. 6. — Results for various ID flow models, showing the log of scaled acceleration w' (left) 
and the ratio of flow speed to escape speed v/v esc (right), plotted vs. scaled inverse radius 
x = 1 — i?*/r. The dotted and dot-dashed curves represent respectively the point-star and 
finite-disk solutions for spherical expansion, while the dashed and solid curves respresent the 
corresponding non-spherical solutions. The panels a-d represent various models, with key 
parameters as labeled on the left. The top three models (a-c) use the Kopp-Holzer scaling 
in eqn. (17), while the bottom model uses the scaling of eqns. (20) and (22), with r/* = 10. 
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found in M-88's solution with gas pressure included), we see that non-spherical expansion 
with h max = 2 leads to a speed increase of 3.46. 

Table 1 summarizes the speed enhancement factors for this well as for a series 

of other models to be described below. 







a = 


0.5 


a = 


0.6 


a = 


0.7 


Model 




r — > oo 


r = 6i?* 


r — > oo 


r = 6i?* 


r — > oo 


r = 6-R* 


M-88: h max 


= 1.5 


1.80 


1.94 


2.01 


2.16 


2.39 


2.56 




= 2.0 


2.22 


2.39 


2.62 


2.81 


3.46 


3.69 


M-88 + f d : h max 


= 1.5 


1.26 


1.25 


1.37 


1.36 


1.60 


1.58 




= 2.0 


1.47 


1.46 


1.71 


1.69 


2.22 


2.20 


KH+UO-02: 77* 


= 0.1 


1.10 


1.08 


1.13 


1.12 


1.21 


1.19 




= 1.0 


1.26 


1.23 


1.38 


1.33 


1.65 


1.56 




= 10. 


1.43 


1.36 


1.67 


1.56 


2.22 


2.01 


eqn.(20): 


= 0.1 


1.13 


1.11 


1.19 


1.15 


1.29 


1.24 




= 1.0 


1.26 


1.21 


1.38 


1.30 


1.64 


1.50 




= 10. 


1.45 


1.34 


1.71 


1.51 


2.33 


1.91 



Table 1: Ratio of flow speeds between non-spherical vs. spherical expansion models. The top 
group of M-88 models is based on the CAK point-star model, but all remaining data include 
the finite-disk correction. The third group uses the Kopp-Holzer non-spherical function (17) 
with parameters set by fits from fig. 7 of UO-02 (see also fig. 7 here), while the bottom 
group uses eqn. (20) with R c = -qV % R*. The boldface values emphasize the strong speed 
enhancement of the M-88 point star model, and the more moderate values of the third and 
fourth model group. Note the latter agree quite well with the MHD simulation results from 
fig. lb, which give polar speed ratios of 1.12, 1.34, and 1.54 for 77* = 0.1, 1, and 10. 

A key point here is that the strong flow speed increases M-88 obtained from including 
non-spherical expansion in the point-star CAK model are quite analogous to the speed 
increases obtained by including a finite-disk correction. Let us next extend our analysis to 
include both effects. 
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4.2. Combined Effects of Non-Spherical Expansion and a Finite Disk 

For a radial flow the combined modification to the line force due to finite disk and 
non-spherical expansion is given simply by the product of these separate 4 correction factors, 

9fd,nse = fdh a gcAK ■ (19) 

We can thus analyze the combined dynamical effect by solving the equation of motion given 
again by eqn.(12), applying now the simple substitution / — > (f d / f*)h a . With this combined 
rescaling, fig. 4 again gives the dependence of \fw' on the combined correction factor /. 

The solid curves in fig. 6a give the acceleration and velocity for the full solution of 
combined non-spherical expansion plus finite-disk. Compared to the above point-star case, 
note the acceleration is now very strong, with indeed the velocity in this case going off the 
scale, approaching an extaor dinar ily large terminal speed = 13.6 v esc = 13,255 km/s! 
However, note that even the spherical, finite-disk model in this case also has a very large 
terminal speed of = 6.13 v esc = 5971 km/s. As such, in relative terms of non-spherical 
to spherical enhancement, the finite-disk case implies a terminal speed factor of "only" 2.22, 
which is actually significantly less than the factor 3.46 of the M-88 point-star assumption 
for this case. 

Nonetheless, in absolute terms the implied terminal speeds well in excess of 10,000 km/s 
in such a model are several times higher than fastest speeds (ca. 5000 km/s) inferred from 
blue-edges of the P-Cygni wind lines from "normal" OB stars (Prinja & Howarth 1986; 
Howarth and Prinja 1987). In this sense, if finite-disk effects had been included in M-88's 
original analysis, the implications of assuming hot-star magnetic fields would have seemed 
even more dramatic; indeed, the lack of observational evidence for such high speeds may 
well have been interpreted as strong evidence against existence of a dynamically significant 
magnetic field (i.e. with confinement parameter 77* > 0.1). 

But now, with the benefit of the full MHD simulations of UO-02, we have a strong 
incentive to investigate the situation further, and indeed the simple formalism here makes it 
possible to readily examine the effect of varying parameters, as summarized in table 1. 



4 Actually, a faster-than-radial expansion generally also implies a larger lateral velocity gradient, and thus 
alters the finite-disk correction factor from its usual spherical form (11). In practice, we find this to be a 
minor correction, except very near the stellar surface, where it leads to a slightly increased mass loss rate; 
see Appendix B. 
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4-2.1. Effect of More Modest CAK Exponent a 

The analysis in §3, and particularly the results in figs. 4 and 5, suggest that a key 
parameter in influencing the wind flow speed is the CAK power- index a. Thus, in comparing 
the M-88 results with those of the UO-02 MHD simultions, it is important to note that the 
former assumed a = 0.7 and the latter a = 0.6. 

Fig. 6b thus provides results for the same Kopp-Holzer parameters as used in M-88 
(and in fig. 6a), but now with a = 0.6. Note that the accelerations and speeds are all 
reduced. In particular, the case with non-spherical expansion combined with finite-disk now 
has a terminal speed of = 6.28t> esc = 6113 km/s, representing a more modest factor 
1.71 enhancement over the terminal speed = 3.67v esc = 3570 km/s of the corresponding 
spherical, finite-disk model. 

4-2.2. Using Dynamically Motivated Values for the Kopp-Holzer Parameters 

To further facilitate comparisons with our MHD simulations results, let us next consider 
Kopp-Holzer parameter values that are adjusted to fit the non-spherical expansion found in 
these simulations. Using eqn. (16) to convert magnetic field variations into flow area, 
fig. 7 plots the radial variation of the polar expansion factor found in the MHD simulations 
assuming various magnetic confinement parameters 77*. The dashed curve shows the expected 
result for pure dipole in limit of large confinement parameter, 77* — > 00. The lower, middle 
and upper solid curves show the MHD results for for 77* = 0.1, 1, and 10. The dotted curves 
plot the best-fit Kopp-Holzer parameterizations, quoted in fig. 7 of UO-02. 

For the Kopp-Holzer fit to the strongest confinement case 77* = 10, fig. 6c shows the 
enhancements in acceleration and speed are again quite modest. Indeed, at the radius 
Rmax — 6-R* of the outer boundary of the UO-02 MHD simulation, the speed of this ID 
model is vq = 5.72v esc = 5576 km/s. Relative to the speed of the spherical case at this 
radius, vq = 3.29v esc =3200 km/s, this respresents an enhancement factor of about 1.56, 
which is now remarkably close to the factor ~ 1.54 found for the polar flow of the 2D MHD 
model (cf. fig. lb). 

Table 1 shows similarly good agreement between this modified Kopp-Holzer parameter- 
ization and the MHD simulations (fig. lb) that show polar speed enhancements of ca. 1.12 
and 1.32 for confinement parameters 77* = 0.1 and 1. 



- 19 - 




Fig. 7. — Non-spherical expansion factor h(r) for polar flow from full MHD simulations of 
UO-02, plotted vs. radius for magnetic confinement parameters rj* = 1/10, 1, and 10 (solid 
curves). The dotted curves compare the best-fit Kopp-Holzer parameterizations, quoted in 
fig. 7 of UO-02 as respectively h max = 1.23, 1.73, and 2.43 and a/R* = 0.73, 0.89, and 1.13, 
with R\ = R* in all cases. The dot-dash curves here add the simple physical fit form of eqns. 
(20) and (22). The dashed curve shows the expected result for pure dipole in limit of large 
confinement parameter, 77* — > 00. 
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4.3. A Dynamically Motivated Non-Spherical Expansion Function 

The overall variation of the MHD non-spherical expansion factors plotted in fig. 7 (and 
in fig. 7 of UO-02) suggests a simpler, more physically motivated fitting-function, 

^> = £tt£- (20) 

which has just a single parameter, R c , representing a "magnetic confinement radius". This 
characterizes the transition between the dipole scaling h ~ r that applies for a strong field 
near the surface, and the asymptotic monopole scaling h — > constant that applies at large 
radii, where the wind has stretched the field into a nearly radial configuration. Note that in 
the opposite limits R c — > or R c — > oo the monopole or dipole scalings apply throughout the 
flow, as appropriate to the limits of weak (77* — > 0) or strong (i]* — > 00) magnetic confinement. 

This suggests a direct scaling of R c with the confinement parameter i]*. In general, 
the ratio of magnetic to wind-kinetic energy should scale as i](r) = 77* (r / R*) 2 ~ 2p , where 
B{r) ~ l/r p , with p = 3 for a dipole, and p = 2 for a monopole. If we then assume the 
confinement radius occurs where this ratio has some order-unity value i] c , we obtain a general 
scaling of the form 



R 
R 



C 



Vc 



(21) 



where q = l/(2p — 2) should be intermediate between the dipole (p = 3) value qd = 1/4 
and the monopole (p = 2) value q m = 1/2. Taking the simple case that the order-unity 
constant r\ c is in fact unity, we find quite good fits for an exponent q = 3/8, representing the 
"geometric mean" between the dipole and monopole scalings for R c , 

R c w vl /8 R* ■ (22) 

The dot-dash curves in fig. 7 plot the radial variations of this simple, physically motivated 
form of the single-parameter non-spherical-expansion function (20). Note that the agreement 
with the corresponding MHD simulation results (solid) is almost as good as the best-fit, 
three-parameter, Kopp-Holzer form [dotted; eqn. (17)] that has less physical justification. 

Fig. 6c shows the acceleration and speed when this simple scaling is applied to the strong 
confinement case with R c = 10 3//8 i?* = 2.37-R*. At the simulation outer radius R max = 6R*, 
the speed ehancement of this non-spherical, finite-disk model is now a factor of about 1.51, 
which is again remarkably close to the factor ~ 1.54 found for the polar flow of the 2D 
MHD model (cf. fig. lb). Table 1 shows similarly good agreement between this simpler 
non-spherical enhancement factor and the MHD polar speed enhancements for the other 
confinement parameter cases. 
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Mass Loss for Tilted Flow 



In addition to modifying the flow divergence, a magnetic field can also alter a flow's 
direction. In particular, for a strong field the initial flow off the stellar surface occurs at an 
angle tilted to the local radial direction. In this section we examine how such a base flow 
tilt should influence the local mass flux. 

Since the mass loss is set very near the stellar surface, where combined factor from non- 
spherical divergence and finite-disk correction has its minimum, let us ignore both sphericity 
and magnetic divergence effects and consider the field and flow to be along a fixed direction 
s, assumed to be tilted by an angle 9 B with respect to the local vertical direction, defined 
by unit vector z. We then have /!# = cos#b = s • z, and in terms of the local horizontal 
direction x, we can also write s = fi B z + a/1 — p 2 B x. 

For our usual simplification of negligible gas pressure, the steady-state equation of mo- 
tion along the fixed field direction s becomes [cf. eqn. (1)] 



(V • V)(§ • V) = -g*(s • Z) + S • gcAK , 



(23) 



where on the right- hand- side gcAK represents the vector radiative force, and we have used 
the fact that surface gravity <?* = (1 — T)GM j R\ is along the local vertical z. On the 
left-hand-side the advective derivative of velocity along s can be written as 



(v-V)(s-v) = v s (s-W)(v s ) 

dv s 



dv s 
dx 



(24) 



where the last expression assumes that there are no horizontal variations in velocity. 



For the line-driving, let us again assume a point-source approximation with the radiation 
confined to a pencil beam along the local vertical z. In terms of an associated surface flux 
F* = L/4-irRl, the projected CAK line- force term then becomes [cf. eqn. (2)] 



s • gCAK = 



1 Q^eF* 



1 — a c 



z • Viz • vl 



I — a c 



1 — a c 
1 — a c 



Qn e pc 
fi B dv s /dz 

Qn e pc 
HBVsdvs/dz 

QK e m s c 
fi B v s dv s /dz 
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where the third equality uses mass conservation along the tilted flow to eliminate the density 
p in favor of the mass flux in the tilted direction, m s = pv s . The final equality then uses 
the fact that the associated radial mass flux component rh r (which sets the associated mass 
loss through the spherical surface) is reduced by another factor p B , i.e. fn r . = pu B m s . After 
dividing through by /is and making the modified definition [cf. eqn.(3)], 

. v s dv s /dz 

w = 

9* 

the momentum equation (23) can again be recast in the familiar form [cf. eqn. (4)], 

w' = -l + Cfi 2 B a {w') a , (27) 

where C is defined in eqn. (5), and is related to the maximal value of the radial mass flux. 
Equation (27) is mathematically identical to eqn. (4) under the substitution C — > C p 2 B . As 
such, following the same procedure as in §2.1, we find the surface mass flux scales as 

m r = Am CAK (28) 

where tHcak = Mqak/^R 2 is the analogous surface mass flux in the CAK model for 
spherical mass loss. 

The p: B reduction of the radial surface mass flux predicted by this simplified planar, ID 
flow analysis is actually in quite good agreement with results of our full 2D MHD simulations 
(Ud-Doula 2002; UO-02). Indeed, the solid curves of fig. la above show that the mass flux 
scaled by p 2 B is roughly constant in latitude, and equal to the mass flux in the spherical case. 
Physically, one factor of \± B can be attributed simply to geometric projection of the tilted 
flow onto the vertical normal to the star's spherical surface, while the other can be seen to 
result from the tilted flow acceleration having less effect in desaturating the optically thick 
absorption of the radially streaming radiation. 

This lower initial mass flux combined with the faster-than-radial area expansion leads to 
a lower density, which through the CAK line-force leads to a greater acceleration and thus a 
faster terminal speed. In cases with a large magnetic confinement parameter (i.e. 77* = 10), 
the simulations thus tend to show the greatest terminal speeds at latitudes just away from 
the magnetic equator, which originate from the open field lines that are furthest from the 
magnetic pole, and for which the initial flow tilt at the surface wind base is greatest. (See 
fig. lb here and right panel of fig. 8 in UO-02). 



(26) 
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6. Discussion 

The mass flux and speed variations arising in such magnetic models of line-driven winds 
provide a potentially attractive mechanism for explaining the extensive wind structure and 
variability commonly inferred from observations of UV P-Cygni lines from OB stars (e.g., 
Kaper et al. 1996; Prinja et al. 1998). If the models here were extended to include even 
a modest stellar rotation along an axis that is tilted relative to the magnetic dipole axis, 
then the speed variations with magnetic latitude would naturally lead to Co-rotating Inter- 
action Regions (CIRs), much as occur between faster and slower streams in the solar wind. 
Following initial suggestions by Mullan (1984), such CIRs now represent a popular general 
paradigm for explaining the commonly observed "Discrete Absorption Components" (DACs) 
in P-Cygni absorption troughs. But in previous hydrodynamical simulations by Cranmer 
and Owocki (1996), the speed variations leading to CIR formation had to be induced some- 
what artificially, by assuming a modulation of the base wind driving, as might result from, 
e.g., bright or dark spots on the stellar surface. Since the analyses here (and the UO-02 
MHD simulations) show that even modest magnetic confinement parameters can lead to 
quite substantial latitudinal variations in flow speed, a tilted dipole field seems a much more 
likely mechanism for inducing such wind CIRs. 

Note that such CIRs would also tend to limit the highest speeds occuring from the 
most rapidly diverging flow, simply from the truncation by the interaction with the slower 
speed flow. As such the cases with the highest overall flow speeds would be those with a 
strong, rotation- aligned dipole, for which the high speed polar wind wouldn't be brought 
into interaction with the slower, equatorial wind. Such cases may be relatively rare, but are 
possible. The flow speeds would be quite high, perhaps up to twice that of a spherical wind. 
But the analysis here shows that they needn't be absurdly high, at least for moderate values 
of the CAK power index, i.e. a 0.6. 

7. Concluding Summary 

This paper analyzes the effects of the flow tilt and non-spherical expansion associated 
with magnetic channeling on the mass flux and flow speeds of a line-driven stellar wind. Our 
main results are summarized as follows: 

1. A faster-than-radial expansion leads to an enhancement in wind speed, but the relative 
corrections are typically about 50% when compared to a finite-disk-corrected model 
with a moderate CAK power-index a = 0.6, much less than the factor three or more 
inferred by M-88 in their analysis of a point-star model with a = 0.7. The analysis 
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here is in good quantitative agreement with results from numerical MHD simulations. 

2. The non-radial expansion obtained from our numerical MHD models differs from the 
heuristic form used by Kopp and Holzer (1976), with the rapid expansion beginning 
right at the wind base, and extending over a quite large radial distance. We propose a 
simpler, physically motivated fit function [eqn. (20)] controlled by a single parameter, 
R c , representing a characteristic radius for wind magnetic confinement, and scaling as 
R c ~ rj^ 8 with the magnetic confinement parameter i]* (= B\ q R\j Mvoo). 

3. The radial mass flux at the stellar surface with a tilted magnetic field is reduced by the 
square of the cosine of the tilt-angle (fi 2 B ) compared to a non-magnetic, spherical wind. 
Physically, one factor of fi B stems from the geometric projection of the tilted flow onto 
a surface normal, while the other arises from the reduced desaturation of absorption 
of radially streaming radiation by acceleration along this tilted flow. 

4. A perturbation analysis (in Appendix A and B) shows that the corrections from a 
small, but non-zero gas pressure scale with the ratio of sound-speed to escape-speed, 
a/v esc . Relative to a zero-sound-speed, finite-disk-corrected spherical wind, typical 
increases in the mass loss rate are 10-20%, with comparable relative decreases in the 
wind terminal speed. In non-spherical models, the additional gas pressure correction 
is typically just 1-2%. 

Overall, the results confirm that even modest magnetic fields can have a substantial 
influence in line-driven winds, with faster areal divergence enhancing the wind acceleration 
and flow speed. But even the largest inferred enhancements still generally allow for speeds 
that are within the range of observational constraints. As such, strong magnetic fields in hot- 
star winds should not be precluded, as they might have been if previous analyses implying 
large speed enhancements had not been reduced by the more complete study here. 
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A. Non-Zero-Sound-Speed Correction for a Spherical, Finite-Disk Model 

Let us examine here the effect of a finite gas pressure associated with a non-zero sound 
speed. In the zero-sound-speed limit considered in §§3-5, the minimal value of the finite- 
disk-correction factor, f d * = 1/(1 + a), occuring near the stellar surface, r — > R*, sets the 
critical conditions that determine the maximal allowed steady mass loss rate, Mfd- But 
with a finite sound speed, the reduced effective inertia (1 — s/w)w' in the base region near 
the sonic point shifts this critical location for most difficult driving slightly away from the 
surface wind base. But then, because of the somewhat larger finite-disk factor, this in turn 
allows a somewhat larger maximal mass loss rate. 

To analyze these effects, let us recast the equation of motion in the form, 

F(x, w, w') = (1 - s/w) w' + l- ^—L—^-C c (w'r - 2s/ (1 -x) = 0, (Al) 
where s = (a/v esc ) 2 is typically of order 10~ 3 , and 

5m s ee - 1 (A2) 

is the fractional correction of the mass loss rate due to this small, but non-zero s. The 
new critical point where this mass loss is fixed now occurs somewhat away from the surface, 
x c > 0, but in the usual cases that s <C 1, we expect both 5m s <C 1 and i c <1. 

In the narrow region from the sonic point to this critical point, the driving is dominated 
by the need to overcome gravity. This implies an outward acceleration that scales as the 
inward gravity, namely w'(x) ~ w' c — a/(l — a), which in turn implies that the velocity 
in this initial region nearly follows the standard CAK ("beta=l/2") solution, w(x) xs w' c x. 
Using this, we obtain a spatially explicit form for the finite-disk factor, which upon expansion 
to first order in x < x c -C 1 becomes 

f(x) « 1 + 4x . (A3) 
We then find the equation of motion (Al) takes the spatially explicit form 

F(x) = ± + Ax-a5m s + 2s = 0, (A4) 

x 

where we have also expanded 1/(1 + 5m s ) a and kept terms to first order in the smallness 
parameter s. The critical point is located where this function has a spatial minimum, 







dF 

dx 



4-^, (A5) 

x 2 c 
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which solves to 

x c w - ajs/2 . (A6) 

Applying this in eqn. (A4) and keeping only the leading term of order y/s, we find the mass 
loss correction is given by 

5m s ~ 4v/1 a — . (A7) 

^ ^esc 

Beyond this critical point, the increase in the finite-disk factor means that solutions to 
the equation of motion quickly approach their large gradient limit, w a ~ a fC c /(l + 5m s ) a . 
(See dashed curves in fig. 4.) Since the wind terminal speed scales as ~ \/v7, this implies 
that the fractional speed corrections from a finite sound speed should scale as 

*W = - 1 « (1 + Sm.)-*™-* - 1 « « -^L= — . (A8) 

^00,5=0 2(1 - a) VI - a f esc 

For a typical case with a — 1/2 and s = (a/v esc ) 2 = 0.001, we find x c 0.012, 5m s ps 0.18, 
and 5t>oo,s ~ —0.09. For the same s but ct = 0.6, we get 5m s 0.13 and <k>oo iS —0.10, 
while for a = 0.7, we find Sm s rs 0.10 and <5foo iS rs —0.12. Applying these corrections to the 
scalings from the zero-sound-speed calculations in §3 generally gives very good agreement 
(to about a percent) with the full solutions. 

It is worth noting here that, even though the corrections for finite gas pressure are 
modest, they are not, perhaps, as small as one might have expected. That is, even though the 
ratio of internal energy to escape energy is of order s = (a/v esc ) 2 , the corrections scale with 
y/s = a/v esc . Moreover, this scaling is multiplied by a large coefficient, which stems from the 
relatively steep increase in the leveraged finite disk factor - e.g. scaling as f 1 ^ a {x) m 1 + 8rc 
for a = 1/2, equivalent to an r 8 dependence on the radius r. As a result, even though the 
reduced inertia of a mildly supersonic flow only leads to a relatively modest outward shift of 
the critical point from the stellar surface, the stronger finite-disk factor leads to a stronger 
radiative driving, and thus allows a larger mass loss. Finally, that larger mass loss also 
implies a comparable relative reduction in the wind terminal speed. 
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B. Non-Zero-Sound-Speed Correction for a Non-Spherical, Finite-Disk Model 

In the case of non-spherical expansion, the analysis follows much as above, except that 
now we must include the non-spherical expansion factor h(x). In the region near and below 
the critical point i < i c C 1 this expands to 

h(x)pal + ti x, (Bl) 

where, for the form defined in eqn. (20), 

K = 1 ^W*ll{l+r 1 fl*). (B2) 

This non-spherical expansion implies a stronger lateral velocity gradient, and for our assumed 
velocity law scaling (w ~ x) this can be shown to modify the finite-disk correction expansion 
(A3) near the surface to 

f h (x) » l + (A + h' )x. (B3) 

Including then also the usual effect of this non-spherical expansion on the density, we make 
the replacement / — > fjji a ~ 1 + [4 + (1 + a)h' ]x and thereby find that equations (A4) and 
(A5) still apply if we simply make the substitution 4 — > 4 + (1 + a)h' Q . The critical point 
thus now occurs at 



'1 — a)s 
4 + (1 + a)h' Q 

while the mass loss correction is now given by 

2 

a 

and the speed correction by 



xc<*\i - A r„ (B4) 



5m s « - y/(l - a) [4 + (1 + a)h' ] — , (B5) 



Svoo s « . . (B6) 

' v/(l-a)[4+(l + aK] 

For the weak, moderate, and strong magnetic confinement cases (??* = 0.1, 1, and 10) we 
find h' ~ 0.3, 0.5, 0.7. For the typical s of order 10~ 3 , we find the extra changes from 
non-spherical divergence range from one to two percent in both the mass loss rate and flow 
speed. In particular, applying the parameters s = 0.0014 and a = 0.6 used in the UO- 
02 MHD simulations, we estimate extra mass loss enhancements for the weak, moderate, 
and strong-field cases to be respectively 0.9%, 1.5%, and 2.1%. These values are roughly 
consistent with the polar mass fluxes plotted in fig. la. 
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